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Abstract 

In this paper, we derive simple and general perturbative formulas for the flavor flux ratios 
Rap = ^va/^vn that could be measured at neutrino telescopes. We discuss in detail the 
role of the uncertainties of the neutrino mixing parameters showing that they have to 
be seriously taken into account in any realistic discussion about flavor measurements at 
neutrino telescopes. In addition, we analyze the impact of such uncertainties in telling the 
standard neutrino oscillation framework from scenarios involving, e.g., neutrino decay and 
we find that the ratio R e ^ is the most sensitive one to "new physics" effects beyond the 
Standard Model. We also compute the more realistic muon-to-shower ratio for a particular 
configuration of the IceCube experiment, observing that using this experimental quantity a 
clear separation between standard and non-standard neutrino physics cannot be obtained. 

1 Introduction 

The standard framework of neutrino oscillations is a successful description of data from 
recent neutrino experiments. However, there might be subleading effects that are not cov- 
ered by neutrino oscillations, since the experimental data are still impaired by rather large 
uncertainties. Such subleading effects could be described by introducing an extended neu- 
trino oscillation framework including so-called damping effects that could stem from, e.g., 
neutrino decay or neutrino decoherence. In this work, we will introduce damping factors, 
describing the damping effects, in a phenomenological way as additional factors to the ordi- 
nary terms in the formulas for the neutrino oscillation probabilities. A promising situation 
to look for such subleading effects would be from neutrino sources at large distances, since 
in this case the result of the damping factors will be most visible because of averaging. 
In addition, it could be an alternative way to measure the fundamental neutrino parame- 
ters, but not the neutrino mass-squared difference, since the averaged neutrino oscillation 
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probabilities do not depend on them. Therefore, we will investigate neutrino oscillations 
with damping effects in general, and different scenarios in which damping effects arise 
from neutrino decay in particular. We will study these scenarios in the case of neutrinos 
coming from large distances. The most plausible way to measure such neutrinos would 
be at neutrino telescopes in ice [1] or water [2]. Many papers in the literature have been 
devoted to study the dependence of the neutrino fluxes on the neutrino mixing parame- 
ters, namely the three leptonic mixing angles #12, #13, and #23 as well as the CP- violating 
phase 5 (for an incomplete list of references, see Refs. [3-9]). The usual starting point is 
that neutrinos are produced via decays of pions and kaons created by hadronic (i.e., pp 
collisions) and photohadronic (jry) interactions [10], leading to the well-known flux ratios 
at the source <^ e : <$P V : </>° t = 1:2:0. However, other effects, like muon energy loss in 
strong magnetic fields [11, 12] can noticeably alter the flavor composition; in addition, a 
different flux ratio at the source can be obtained for the so-called "neutron beam sources" 
in which neutrinos are produced from neutron decays [13,14]. In any case, on their way 
from the source to the Earth, neutrinos oscillate and the fluxes arriving at the detector 
acquire a dependence on the neutrino mixing parameters which can, in principle, be used 
to improve our knowledge on the fundamental neutrino parameters and/or used in connec- 
tion with reactor experiments and neutrino beams to better the determination of S and the 
neutrino mass hierarchy. However, it could happen at the time when neutrino telescopes 
will be operational that the uncertainties of the neutrino mixing parameters could be large 
enough not to be neglected in theoretical estimates of the arriving fluxes [15]. This implies 
that any possible dependence on the neutrino mixing parameters can be completely hid- 
den, and what is also important, that any small deviation from the expected number of 
neutrinos can simply be the effect of our ignorance about the precise values of the neutrino 
mixing parameters and not due to any "new physics" effects in large distance neutrino 
oscillations, as those given, e.g., by neutrino decay [16, 17], breakdown of fundamental 
symmetries [14,18-20], or pseudo-Dirac nature of neutrinos [21,22]. 

For this reason, in this work, we analyze in some detail the impact of the current 
neutrino mixing parameter errors on the determination of the fluxes arriving at neutrino 
telescopes. We explicitly illustrate the dependence on the deviation from the best-fit values 
of 0i2,0i3, and 623 up to second order, in the case that the fluxes at the source are in the 
ratios <^ e : 0° : <f>® T = 1:2:0. For comparison, we also show how the fluxes on Earth 
change if neutrino decay takes place and we discuss whether it is possible to distinguish 
them from the fluxes computed in the standard neutrino oscillation framework, and to 
what extent. 

Our paper is organized as follows: In Sec. O we derive the averaged neutrino oscillation 
probabilities including damping factors. Next, in Sec. El we analyze the three-flavor flux 
ratios and study the impact of the uncertainties of the fundamental neutrino parameters. 
Then, in Sec. IH we investigate the flux ratios with neutrino decay, present the different 
neutrino decay scenarios, as well as we discuss the possibility to detect the subleading 
effects at neutrino telescopes. Finally, in Sec. EJ we summarize our work as well as we 
present our conclusions. 
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2 Derivation of averaged neutrino oscillation prob- 
abilities including damping factors 

In general, the neutrino oscillation probabilities in the standard three-flavor framework 
have the following form: 

3 3 

i=l i=j 

where J^g = UajU^-U^Upi and $y = Am|l/(2£). Here U is the leptonic mixing matrix, 
Am| = m? — m 2 j are the neutrino mass-squared differences, L is the baseline length, and 
E is the neutrino energy. The standard parameterization of the leptonic mixing matrix is 
given by 

C12C13 S12C13 s 13 e~ lScpS 

U = | -S12C23 - Ci2Si 3 S23e 1<5cp Ci 2 C 2 3 ~ Si 2 Sl3S23e 1<5cP Ci 3 S 2 3 J ■ (2) 

•S12S23 - Cl 2 Sl3C23e 1<5cP -Cl 2 S 2 3 ~ Sl2Sl3C23e 1<5cP C13C23 

where Cy = cos 6ij, Sij = sin^, and 5cp is the Dirac CP-violating phase. In the case 
that the ratio L/E is large, i.e., L/E S> 2(Amf •) , one obtains in vacuu the following 
averaged neutrino oscillation probabilities 



{p aP ) = Y. J % = Y,\ u ^\ 2 \ u ^ 2 - ( 3 ) 



Note that these probabilities depend only on the parameters Jjg, where i = j, as well as 
they are independent of the oscillation frequencies, since these have been averaged out to 
zero for large values of the ratio L/E. 

Now, introducing damping factors in Eq. ((T]), the neutrino oscillation probabilities 
can be written as 

3 3 

i=l j=l 
3 

= £A4 + 2 E A.lJ^lcos^ + argJ^), (4) 
i=l 1<«J<3 

where the damping factors are given by 

/ |Am?,|^\ 
Dij = exp I -a>ij — I (5) 

with ctij being elements in a non-negative damping coefficient matrix. We will assume 
OLij to be independent of energy. Note that both the oscillation frequencies and the 
damping factors Dy are functions of the baseline length L and the neutrino energy E. The 
oscillation frequencies are functions of the ratio L/E only, whereas the damping factors 
are more general functions of the parameters L and E, i.e., = Dij(L@ /E 1 ). 

Similarly, as in the averaging procedure above for the standard framework, one can 
perform an averaging of the neutrino oscillation probabilities including damping factors 
for large values of the ratio L/E, i.e., L/E S> 2(Am-j) . Note that it would be unnatural 



1 Matter effects inside the source can affect the neutrino transition probabilities, see Ref. [23]. 
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to include damping effects after one has performed averaging. Thus, in the case when 
(3 = 7, the damping factors are functions of (L/E) 13 only, and therefore, the following 
averaging (with i = L/E) can be introduced^] 

{Pap) = lim f ° XP ^ de = lim i f P^(/) d£ (6) 
x ~^°° Jo ®" x ~^°° x Jo 

In this work, we will not consider the case when (3 ^ 7. Partly, because this case has less 
obvious applications. In general, for an arbitrary value of the parameter (3, it is not possible 
to compute the averaged neutrino oscillation probabilities including damping factors, but 
it is possible for the cases when (3 = 1 and (3 = 2, which are the cases that are important 
for applications to scenarios that could arise in Nature, see Ref. [24]. 
In the case when /3 = 1, we obtain 

3 3 

00 1=1 j=i 

where the effective damping factors are given by 



= -» ( g ) 



1 — exp ( — i x _ aij\AmfA^x 



which can be written in a more simpler form as 

1 (-aijx) 1 



WW-E^tt^. (•) 

n=0 

where 

Am?. 

(Hj = i— 2^ + aii|Am^| ? . (10) 

For £ 7^ 0, the parameters a^- are in general non-zero if i 7^ j, and therefore, one finds that 
\\m. x ^ oa {Dij){x) = 0, whereas a^- = if i = j, which means that (Da){x) = 1, and hence, 

one regains the formula in Eq. ©. On the other hand, for £ = 0, a^- = a^- + i — ^ if * 7^ J , 
and especially, a^- = ajj if i = j. Thus, in the case that £ = and in the limit x — > 00, we 
find that {P a p) = 0, since linx r ^ 00 (.Dj ? )(x) = for all x if % = j and ]im a ._ >0C (.Dy)(a:) = if 
i 7^ j and «jj 7^ 0. However, if = for fixed i = j, then lim I .^ 00 (Z)jj)(x) = 1. Thus, in 
the case that £ = and an = for fixed i, the averaged neutrino oscillation probabilities 
become 

(Pap) = J% = \U ai \ 2 \U pi \ 2 . (11) 
In general, if an = for more than one fixed i, then one obtains 

(P^) = ^J^ = ^|c/ m | 2 |^| 2 , (12) 

where N can be any of the sets {1,2}, {1,3}, {2,3}, and {1,2,3}. In the case when 
N = {1, 2, 3}, one recovers the standard framework formula in Eq. ([3"j). 



2 Note that averaging with respect to large distances can be obtained by simply assuming E to be a constant, 
i.e., by defining I = -gZ, where E is a constant. 
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Similarly, in the case when j3 = 2, we can again write the formula for the averaged 
neutrino oscillation probabilities including damping factors as in Eq. ([3), but now the 
effective damping factors are given by 



(Dij){x) 



1 



2 Jc^-j x 



7r exp 



(Am 



( & m ij' 
erf I kijX + i- 



4k 



ierfi 



v 



4ka 



(13) 



where k 



ij 



aij\Am?j\£ . Note that the argument of the imaginary error function is 



independent of x. Again, the situation is the same as in the case when ,9 = 1, which means 
that non-zero averaged neutrino oscillation probabilities can only be obtained for £ = 
and ouj = for fixed i. 



3 Neutrino flux ratios 

In the cases /3 = 7 = 1 and j3 = 7 = 2, the averaged neutrino oscillation probabilities with 
damping factors included that are described in Eq. (fT2l) can effectively be written as 

3 

(Pap)=J2diJ^, (14) 

8=1 

where di are the normalized damping factors, which can have the value or 1. Of course, 
if di = 1 for i = 1, 2, 3, then one finds the averaged neutrino oscillation probabilities in the 
standard framework, which are given in Eq. ©. 

Starting from the flux ratios at the source (j)® e : <$P V : (\P Vt = 1:2:0, the neutrino 
fluxes arriving at a detector are sensitive to neutrino oscillation in vacuum and are then 
computed as 

<t>v e = (Pee)+2{P^ e ), 

K = <^)+2(P MM >, (15) 
<Pu T = (P er > + 2(iV>. 
Then, the three-flavor flux ratios analyzed in this work are defined as follows: 

9^ 9ur 9u T 

Although the flux ratios are not independent of each other, since a measurement of two of 
them will give the value of the third, we prefer to discuss them separately. 

3.1 Standard neutrino framework 

In this subsection, we analyze the analytical form of three-flavor flux ratios on Earth. These 
ratios depend on three mixing angles and one CP-violating phase, as implied by Eq. (|15l) . 
and are, in principle, sensitive to subleading effects (encoded in the dependence of #13) and 
the errors of the mixing angles 612 and 623- For maximal mixing of #23, i.e., 623 = vr/4, 
and vanishing #13, the robust prediction of Eq. (fT5|) is (j) Ue : (f> v : <f) VT = 1 : 1 : 1 [3], 
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independently of 0i2- Going beyond this zero-order approximation allows us not only to 
study the role of non- maximal 023 and non- vanishing #13 , but also to understand the impact 
of introducing the uncertainties of the neutrino mixing parameters. 

Since the exact formulas for the flux ratios are quite cumbersome, we prefer to present 
our results in some useful approximations. From experimental data, we know that #13 is 
small and 623 is close to maximal mixing. For that reason, we can expand the ratios in 
#13 and 823 = 023 — 7r /4, parametrizing the deviation from maximal 23-mixing [25]. At the 
same time, we can also expand for small 612 = 612 — #12, #12 being the best-fit value for 
9\2- No restrictions have been applied on the CP-violating phase 5, which means that the 
following formulas are valid to all orders in 5. Up to second order in the small quantities 
#13, ^12, and 623, they read 

R efM = 1 + ~ eos(<5) sin (40 12 ) 13 - ~ sin 2 (20 12 ) 5 23 



1 



+ - cos 2 (5) [3 cos (40i 2 ) - 5] sin 2 (20 i2 ) 0f 3 
8 

+ ^ [-28 cos (4#i 2 ) + 3 cos (80i 2 ) - 103] 8 2 2Z 
+ 3 cos(<5) cos (40 i2 ) <5i2#i3 - 3 sin (40 i2 ) ^12^23 

+ 4 cos(<5) [3 sin (80iz) - 22 sin (40 12 )] #13<5 23 , (17) 
lb 

Rer = 1 + \ COS(J) Sin (40 12 ) #13 - \ sin 2 (20 12 ) 5 23 

+ \ cos 2 {5) [3 cos (40i 2 ) + 11] sin 2 (20 i2 ) 2 3 

o 

+ i [4 cos (40 i2 ) + 3 cos (80i 2 ) + 121] 5 2 23 



+ 3 cos(d) cos (40 i2 ) <5i 2 0i3 - 3 sin (40 i2 ) 5 l2 5: 



■23 



+ 4 cos(5) [10 sin (40i 2 ) + 3 sin (80 i2 )] 0i 3 <5 23 , (18) 
lb 

RflT = l 

+ 2 cos 2 (5) sin 2 (20 12 ) 2 3 + [cos (40 12 ) + 7] 8% 3 

+ 2 cos((5) sin (40i 2 ) 013^23 • (19) 

Note that another series expansion up to first order with a different parametrization of the 
initial neutrino flux ratios for the neutrino oscillation probabilities has been presented in 
Ref. [9,2b]. Now, some comments of the above formulas are in order. First, we discuss the 
flux ratio R efjL . In Eq. (|17p . we have clearly separated first- and second-order terms in the 
series expansion. The approximate relation is obviously much more accurate if referred 
to small values of 0i3, 5i2, and 823 ■ In order to check how good our series expansion 
is, we show in Fig. [TJ the ratio between the approximate formulas (at first and second 
order, i?g^ st and i?|^ cond , respectively) and the exact formula for R efl7 as a function of 
the small and unknown mixing angle 0i3. Furthermore, note that the accidental sum-rule 
R efl — R eT + R^ T = 1 holds up to second order in perturbation theory. 

In order to fix some values for the fundamental neutrino parameters, we observe that 
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Figure 1: Comparison between the approximate formulas R^ and R s e e ^ ond as well as the exact 
formula for R efl . The solid curve corresponds to the reference value 1, whereas the dotted and 
dashed curves correspond to the first- and second-order approximations, respectively. 



the following 99 % C.L. limits hold [27] 

30° < 0ia < 38° 

36° < 023 < 54° (20) 
< 013 < 10° 

with best-fit values corresponding to 9\2 = 33°, 023 = 45°, and 0i3 = 0. Thus, in order 
to be sensitive to 6\2 and £23, we choose 0i2 = 38° and 023 = 50° (corresponding to 
612 = #23 = 5°). 

It can be clearly seen that the second-order approximation reproduces the exact values 
at the level of 1 % to 3 %, the worst case being obtained for larger values of the mixing angle 
013. The larger discrepancy between the first-order approximation and the reference value 
is to be mainly ascribed to both the relatively large ^23, which means that second-order 
terms in this quantity cannot be safely neglected, and the absence of terms of 0(612) ■ 

At first order in the small parameters, one important feature of this ratio is that it is 
independent of the error of 0i2, which only enter in higher-order terms coupled with the 
error of the mixing angle 023- This means that, once you have fixed the value of the mixing 
angle 0i2, its uncertainties are not relevant for the determination of R efJi . Thus, this ratio 
mainly depends on 0i3 and the uncertainty of 023- 

In order to estimate which of these parameters encodes the largest part of the uncer- 
tainty of Refi, we can evaluate the maximum spread due to the combination x = cos(6) 0i3 
(x = 623) at fixed #23 (cos(<5) 013), namely 

ao __ Re^(Xmax) ~ Re^( x mm) x 

R efi (x = 0) ' 1 ' 

^max.min being the maximum and minimum value of the variable x. It is simple to derive 
Ai? eA1 in the approximations used above. For x = cos(£)0i3, x m ax,min = ±7r/18 [see 
Eq. (I20p ]. and for vanishing 612 and #23, we simply obtain 

AR e ^ = ^ sin 2 (40 12 ) ~ 20 %, (22) 
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Figure 2: The spread of the three flux ratios considered in this work due to the present uncertain- 
ties of the neutrino mixing parameters [see Eq. (G2|)/ ; around the best-fit values 612 = 33° and 
623 = 45° , as a function of the mixing angle #13. From above to below: R eT (red), i? MT (blue), 
and R efJL (grey). 



for 612 at the best-fit value. On the other hand, for x = 623, ^max.min = ±97r/180, we 
obtain 

Ai? eM = — sin 2 (20 12 ) ~ 45 %, (23) 

which means that the error of 623 contains the largest part of the uncertainty of R e ^ [28]. 

Then, we discuss the flux ratio R eT . Due to the approximate fi— r symmetry [29-33] (for 
an incomplete list of references, see Ref. [9]), the flux ratio R eT shares similar characteristics 
as R e ^, as it can be seen in Eq. (|18p . In particular, the first-order approximations are 
exactly the same, whereas differences appear in terms of 0(8f 3 ), 0(5| 3 ), and 0(613623)- 

Finally, we discuss the flux ratio R^ T . This flux ratio is very peculiar. At first order 
in small parameters, it is exactly 1 and the deviation from 1 is only present when second- 
order terms are included. However, no dependence on the uncertainty 5\2 is present at 
least up to 0(^12) 5 which means that R^ T is mainly sensitive to the central value of the 
mixing angle #12, but, once fixed, the uncertainty of it does not affect the flux ratio. What 
is also an important feature, R^r is always larger than 1, as it can be seen considering that 
the terms in Eq. (|19p under the conditions in Eq. (|20p cannot give negative corrections 
to 1. This feature is not lost in the exact evaluation, as we have numerically checked. In 
summary, the flux ratio R^ T is very close to 1, even if the uncertainties of the neutrino 
mixing parameters are taken into account, since the deviation from its standard value 
is an effect of second order in small quantities. For that reason, if some "new physics" 
mechanism would be able to produce a large deviation from 1, the R^ T flux ratio could 
give us an interesting possibility to study such a "new physics" scenario. 

In order to summarize the previous discussion about uncertainties, we have evaluated 
the whole spread of the flux ratios R e ^, R eT , and R^ T due to the present lack of knowledge 
of the neutrino mixing parameters. We have computed the minimum and maximum values 
of each R a p varying 6*12, 6*23, and 5 in the 99 % C.L. of Eq. (|20l) and presented the results 
in Fig. [2] as a function of #13. From above to below, colored areas contain the possible 
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Figure 3: The parameter space giving a value R efl = 1 ± 5 %. The dotted, solid, and dashed 
curves correspond to 9 2 3 = 43°, 45°, 47°, respectively. The mixing angle Q\i is fixed to 33°. 

values for R er (red), R^ T (blue), and R etl (grey), respectively. As it can be easily seen, 
i? e/i can be as large (small) as 1.18 (0.72), whereas R eT can assume values ranging from 
0.83 to 1.40. On the other hand, i?^ T can deviate from unity by 25 %. 

Thus, it seems to be clear that even in the case that the flux ratios could be measured 
with infinite precision at neutrino telescopes, it will be very difficult to estimate or put 
any reasonable bounds on the neutrino mixingangles (and the CP-violating phaseH) in 
the standard three-flavor oscillation frameworko In particular, the ratios do not depend 
in a dramatic way on #13 (both Figs. Q] and [2] numerically confirm this statement) and 
the dependence on 5 is completely overshaded by the uncertainty of the mixing angle #23- 
However, assuming that some long-baseline neutrino experiments gave us a very precise 
measurement of the mixing angles #12 and 623 (which will certainly not happen in ten 
years from now) and assuming also the quite unrealistic scenario in which the flux ratio 
R e ^ will be measured within 5 % precision at future neutrino telescopes, #13 and 5 will 
not be strongly constrained. In Fig. we quantify this statement showing the allowed 
pairs (9i3,5), which give R efl = 1 ± 5 %, for different values of 623, namely 43°, 45°, and 
47°, represented by dotted, solid, and dashed curves, respectively. The admitted parameter 
space is the region on the left-hand side of the curves. In principle, some values of #13 and 5 
can be excluded, depending on the assumed value of the mixing angle 623, but the precision 
obtained cannot be comparable with that usually claimed for neutrino factories and f3- 
beams. However, as already discussed in Ref. [15], the combination of the information 
coming from neutrino telescopes and those from future neutrino experiment could help in 
constraining some of the neutrino mixing parameters, but only in the optimistic hypothesis 
that one can clearly separate the different neutrino sources (i.e., if the flux of high-energy 
neutrinos is under control). 

3 In Ref. [4], the possibility of measuring the CP-violating phase 6 using neutrinos from far distance was 
discussed. 

4 Note that a first criticism about the possibility to measure #13 and S at neutrino telescopes can be found in 
Ref. [34]. 
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4 Flux ratios including neutrino decay 



Neutrino decay has been invoked in the past to solve both the solar and atmospheric 
neutrino data problems (see e.g. Refs. [35-37]). Decay will deplete the flux of neutrinos on 
Earth by the exponential factor 

«•>(-£) — (-it) 

in which r is the rest-frame neutrino life-time, L the distance between the source and the 
detector, and E and m are the neutrino energy and mass, respectively. Since neutrino 
decay scenarios have not been validated so far [38,39], we can roughly estimate the order 
of magnitude of the allowed values of the ratio r jm requiring that |f < 1. Thus, giving 
a lower bound r/m > 10 -4 s/eV, obtained for a typical L/E solar neutrino [40]. Note that 
from the supernova SN1987A neutrino events, one would obtain a stronger lower bound at 
the level of r/m > 0(1O 5 ) s/eV. However, this bound is not reliable if, as it is the case, the 
value of 9 12 is large [41,42]. Given the bounds previously discussed, we cannot eliminate 
the possibility of astrophysical neutrino decay. 

Two-body modes of Majoron type [43, 44] are still viable neutrino decay processes, 
since they are not strongly constrained (contrary, for example, to radiative decay modes or 
three-body modes, see e.g. Refs. [16,45] and references therein). Stringent bounds on this 
class of models have been derived in Ref. [46] , although the robustness of their conclusion 
has been questioned in Ref. [47]. Thus, we retain the possibility for neutrino decay in 
Majoron models to be a possible mechanism for flux ratio modifications. We adopt here 
the simplifying assumptions that the neutrinos completely decay into the lightest mass 
eigenstate or the lightest and next-to- lightest mass eigenstates over very large L/E and 
that there are no detectable decay products. Two different neutrino decay scenarios (with 
three active neutrino flavors) will be investigated. In the first scenario, which will be 
called "neutrino decay I" , the two heaviest neutrino mass eigenstates can decay and only 
the lightest neutrino mass eigenstate is stable, whereas in the second scenario, which will 
be called "neutrino decay II", only the heaviest neutrino mass eigenstate can decay and 
the two lightest neutrino mass eigenstates are stable. The first neutrino decay scenario, 
neutrino decay I, has been heavily discussed in the literature, see e.g. Ref. [48]. Note 
that both neutrino decay scenarios can be considered for both normal (NH) and inverted 
neutrino mass hierarchy (IH), i.e., for the cases, mi < ni2 < 7773, where mi is the mass 
of the lightest neutrino mass eigenstate, and 771,3 < m i < m 2, where 7773 is the lightest 
neutrino mass eigenstate. The different neutrino decays scenarios are schematically shown 
in Fig. HI 

In the case of neutrino decay scenarios, the parameter (3 = 1, and therefore, we can 
use Eq. (I14p in order to compute averaged neutrino oscillation probabilities (including 
damping factors). Note that in the case of neutrino decay, the damping coefficient matrix 
is described by neutrino decay parameters [24], which are defined as ay = (a, + ctj)/2, 
where = rrtj /Vj is the decay rate with mi being the mass of the ith mass eigenstate and Tj 
is its life-time (in its own rest frame) [48]. Thus, in the case of "neutrino decay I", we have 
an = or 033 = and all other ay's are non-zero in general, since only ai = (NH) or 
03 = (IH), which means that d\ = 1 or d% = 1, i.e., only one of the di's is equal to 1, and 
the others are zero (or, in other words, that the exponential factor in Eq. (I24D is unity for 
the lightest mass eigenstate and vanishing for the others). In the case of "neutrino decay 
II", we have an = 022 = or an = 033 = 0, since ai = 02 = (NH) or a\ = 03 = (IH), 
which means that two of the d^s are different from zero, i.e., d\ = c?2 = 1 or d\ = d% = 1. 
It is now easy to obtain flux ratios on Earth for both scenarios. 
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Figure 4: The two different neutrino decay scenarios. Upper-left picture: Neutrino decay I with 
normal hierarchy {an — =>- d\ — 1). Upper-right picture: Neutrino decay I with inverted 
hierarchy {a 33 = =>- d 3 = 1). Lower-left picture: Neutrino decay II with normal hierarchy 
(a n = a 2 2 = =>■ &\ — &i — 1). Lower-right picture: Neutrino decay II with inverted hierarchy 
(an = a 33 = ^ di = d 3 = 1). 

4.1 Neutrino decay I 

The transition probabilities needed to build the flux ratios for the "neutrino decay I" 
scenario, illustrated in the first row of Fig. [4J are obtained from Eq. (|14p in which the 
terms diJ l \ have to be evaluated in the unique index corresponding to the stable mass 
eigenstate (i = 1 for NH and i = 3 for IH). 

For the sake of simplicity and with the aim of checking that our damping formalism 
reproduces the results obtained with the standard formulation of large L/E behavior of 
the flux ratios, we report here only the first-order approximations H Note that the flux 
ratios are independent of the initial flavor composition at the source, and according to 
Ref. [17], we find that <j) Ve : <j> Vfi : <$> Vt = |^ei| 2 : |f^i| 2 : |^ri| 2 for the NH and <f) Ve : cj) Vfl : 
<f>v T = l^esl 2 : l^sl 2 = l^l 2 for the IH. 

First, let us discuss the NH case: 

R efX = 2 cot 2 (0i2 ) 

- 4 cot (6*i 2 ) esc 2 (6*12 ) 8 X2 - 4cos(«5) cot 3 (0i 2 ) 0i3 + 4 cot 2 (0 12 ) 5 23 , (25) 

R eT = 2 COt 2 (6»l 2 ) 

— 4 COt (^12) CSC 2 (012) 8\2 + 4 C0S(<5) COt 3 (012) 013 — 4 cot 2 (0i 2 ) <5 23 , (26) 

Rpr = 1 + 4 cos(<5) cot (0i 2 ) 0i3 - 45 2 3- (27) 

The main feature of these formulas are, except for the flux ratio R^t, that the zeroth-order 
approximation is different from 1, being dependent on cot 2 (#12)- For the best-fit value used 
in this work, this means that <j) Ue : <p Vii : <p Ur ~ 5 : 1 : 1 [17], a very large deviation compared 
to the standard neutrino framework. Moreover, unlike Eqs. (|17p and <\18h . a dependence 
on <5i2 appears at first order, which means that the error of the mixing angle 0i2 cannot be 

5 In Table [TJ the goodness of our approximate formulas can be found. 
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neglected. Note that R^ T obtains corrections with respect to 1 also at first order. Thus, 
we expect that the intrinsic uncertainty of the flux ratios (given by our ignorance on the 
neutrino mixing parameters) is much larger than the standard framework case. 
Second, in the case of IH, the flux ratios assume the following simple structure 

Refj, = Rer = 0, (28) 

R IXT = l + 4<5 23 . (29) 

The signature of this particular case is extremely different from both the standard frame- 
work and the "neutrino decay I" with NH. Corrections to R e ^ = and R er = appear at 
second order in small quantities and are dependent on #13. Thus, we expect a substantial 
deviation from zero at large #13, since corrections appear only at second order in #13 [49]. 

4.2 Neutrino decay II 

Similarly, the "neutrino decay II" scenario is depicted in the second row in Fig. [U In 
this case, the two lightest neutrino mass eigenstates are stable instead of only the lightest 
in the "neutrino decay I" scenario. Therefore, we will have two dtJ™p terms in Eq. (|14p 
instead of only one (i = 1 and i = 2 for NH and i = 1 and i = 3 for IH). In App. IA.H we 
present formulas for the flux ratios only up to first order in the small parameters #13, <5i2 , 
and #23, since the first-order formulas make a much better approximation for "neutrino 
decay II" than "neutrino decay I". Note that for NH, the formulas are independent of 5\2- 
In addition, for both NH and IH, all formulas are dependent on £23, which means that the 
uncertainty of #23 is crucial. Furthermore, we observe that the structure of the formulas for 
NH and IH is very different. For IH, nearly all terms (except from the zeroth-order term of 
the flux ratio R^ T ) include the factor 3 — cos(2#i2) = 2 + 2sin 2 (#i2) in the denominators, 
which has a value between 2 and 4 and is therefore not equal to zero at any time though. 
In addition, R eiL is independent of the CP-violating phase 5, whereas R er is dependent on 
<5, which actually spoils the approximate \i — r symmetry. In some sense, the values of 
the flux ratios are closer to the standard neutrino framework for "neutrino decay II" than 
"neutrino decay I" 

4.3 Discussion about the neutrino decay scenarios 

The previous considerations have been summarized in Figs. [5] and El in which we present 
the flux ratios R efl . and R^ T , respectively, using the two neutrino decay scenarios with both 
NH and BB0. 

In Fig. \E\ for the flux ratio R ef j,, we observe that the "neutrino decay I" scenario for 
both NH and IH can, in principle, be distinguished from the standard neutrino oscillation 
framework. As expected from the very different first-order approximations in Eqs. (|25p and 
(1281) . even including the errors of the neutrino mixing parameters, the bands representing 
the total spread of R efl are very well separated, especially for the IH case in which a 
quadratic behavior of #13 can be clearly seen (notice that we are using a log scale on the 
vertical axis). In addition, we observe the large dependence of #13 on the uncertainty of 
the NH case, which is mainly to be ascribed to the term 4cos(5) cot 3 (#12) &13- 

A very different situation arises for the "neutrino decay II" scenario (solid and dotted 
curves for NH and IH in Fig.[5j respectively). Due to the fact that this scenario contains two 

6 Again, in Table QJ we present the numerical results on the flux ratios. 

7 Note that the flux ratio R eT is very similar to the flux ratio R efJ ,, and we have therefore not presented this 
flux ratio in a separate figure. 
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Figure 5: The flux ratio R eil including neutrino decay as a function of the mixing angle #13. The 
red (middle) area corresponds to the standard neutrino oscillation framework, whereas the blue 
(upper) and green (lower) areas correspond to neutrino decay with the lightest neutrino mass 
eigenstate being stable only for NH and IH, respectively. The case of the two lightest neutrino 
mass eigenstates being stable is visualized by the solid (NH) and dotted (IH) curves, respectively. 
The uncertainties of the neutrino mixing angles are those given in Eq. (EFJj, around the best-fit 
values 612 = 33° and #23 = 45°. 
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Figure 6: The flux ratio R^ T including neutrino decay as a function of the mixing angle #13. 
The red (innermost) area corresponds to the standard neutrino oscillation framework, whereas 
the blue (outermost) and green (middle) areas correspond to neutrino decay with the lightest 
neutrino mass eigenstate being stable only for NH and IH, respectively. The case of the two 
lightest neutrino mass eigenstates being stable is visualized by the solid (NH) and dotted (IH) 
curves, respectively. The uncertainties of the neutrino mixing angles are those given in Eq. $2(h) . 
around the best-fit values Qyi = 33° and #23 = 45°. 

stable mass eigenstates, it is quite similar to the standard neutrino oscillation framework, 
in which all the mass eigenstates are not allowed to decay. Therefore, differences in R efl 
are very difficult to be discerned (the same comment can be applied to Fig. [6]). 

In Fig. [6j the results of the flux ratio R^ T are shown. As already outlined in the 
previous section, in the standard framework, R^ T is always larger than 1, even including 
the uncertainties of the neutrino mixing angles. The neutrino decay scenarios allow i?^ T 
to be smaller than 1 (which could be a clear signature of "new physics" effects), but the 
large spread of the ratio at any value of #13 does not allow an explicit separation between 
the standard and decay phenomenologies. Thus, this flux ratio is not a good candidate to 
find deviations from the standard neutrino oscillation framework. 

In Table [U we present numerical calculations of the flux ratios <p Ue : : <fiu T f° r the 
different neutrino decay scenarios, comparing the zeroth-, first-, and second-order results 
with exact calculations!! In general, we observe that the second-order results contains 
all the features connected with non- maximal #23 an d non- vanishing #13. Thus, it seems 
to be important to include second-order effects in any analytical treatment of flux ratios. 
However, note that the first-order results work especially well for the standard framework 
and "neutrino decay II", whereas for "neutrino decay I", the first-order results do not 
reproduce all features of the exact computations. 

8 Note that the zeroth-order values for the flux ratios in the case of "neutrino decay I" with both NH and IH 
can, e.g., also be found in Ref. [17]. 
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<Pu e ■ 0„„ : <f> Vr 


th 




2 nd 


exact 


standard framework 


1:1:1 


0.93:1:1 


0.93:1.06:1 


0.93:1.06:1 


neutrino decay I (NH) 


4.74:1:1 


3.08:1.06:1 


3.14:0.99:1 


3.20:0.99:1 


neutrino decay I (IH) 


0:1:1 


0:1.35:1 


0.02:1.41:1 


0.02:1.42:1 


neutrino decay II (NH) 


2:1:1 


1.81:0.65:1 


1.78:0.74:1 


1.78:0.73:1 


neutrino decay II (IH) 


1.08:1:1 


0.91:1.28:1 


0.93:1.30:1 


0.93:1.30:1 



Table 1: Flux ratios (j> Ve : (f) Ufi : (p Ur for the different neutrino decay scenarios. , 1 st , and 
2 nd represent the results obtained up to zeroth, first, and second order in small parameters, 
respectively, whereas the exact computation is given in the column exact. The values of the 
fundamental neutrino parameters used are: Qyi — 33° , #13 = 5°, and 5 = 40°. In addition, we 
have assumed 8 12 = 5° and 5 2 3 = 5°. 

4.4 Discussion about neutrino telescopes 

In the previous sections, we have discussed the flux ratios from a theoretical point of view; 
thus taking into account the role of the uncertainties of the neutrino mixing parameters. 
Let us now address the question of what can really be measured at neutrino telescopes 
and how the statistics, i.e., the number of expected events, affects the results that can be 
obtained. 

The measurement of neutrino flavor seems to be a very difficult task, which means that 
the three flux ratios previously studied are unlikely to be directly accessible. However, 
they can be inferred from distinct classes of events: showers, muon tracks, and unique v T 
signatures (see, e.g., Ref. [50]). All neutrino flavors undergo neutral-current interactions 
(inside or nearby the detector), which result in hadronic showers in a detector. On the 
other hand, u e charged-current interactions also produce electromagnetic showers that, in 
principle, could be distinguished from the hadronic ones due to the different muon content 
(which is absent in the electromagnetic showers). The shower rate also includes v T charged- 
current events, since at least for energies below the order of PeV the tau track cannot be 
separated from the shower. For energies above the order of PeV, double-bang [51,52] and 
lollipop events are also possible, but we do not discuss these tau signatures in this work. 
Muon tracks originate from charged-current interactions; muons always emerge from 
the shower, which means that this kind of process can be distinguished from a typical 
shower event. 

Neutrino telescopes seem to be the ideal place where to seek for high-energy neutrino 
interactions. In deep ice or water, neutrinos are detected by observation of the Cherenkov 
light emitted by charged particles produced in charged-current and neutral-current interac- 
tions. Such detectors are mainly sensitive to TeV-PeV neutrino energies; thus opening the 
possibility of studying this still unexplored energy regime. In order to observe a hadronic 
event inside a telescope, a neutrino must first survive as it crosses the atmosphere and 
then the ice (or water) above the detector. Its typical interaction length in a medium of 
density p is Lq = 1/(pNac uN ), where a uN = Gqq + cr^c * s the total cross-section to have 
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an interaction that depletes the neutrino energy. It is usual to express this length in terms 
of its depth: xq = pLo (i.e., one meter of water has a depth of 100 g/cm 2 ). 

A neutrino from a zenith angle 9 Z must cross a column density of material x(6 z ) = 
L dl p(l,9 z ). In practice, the path in the atmosphere is negligible and x(9 z ) is just the 
depth of the ice or water above the detector. The probability that it does not interact 
before reaching the detector is then P S urv(E u ,0 z ) = e ~ x / x o. Once in the detector, the 
probability of an event is P\ n t(E v ) ~ 1 — e~ LpNA<T int , where L is the linear dimension of 
the detector and is the cross-section for the interaction taken into account (charged or 
neutral). Therefore, the number of shower events in the telescope in an observation time 
T for a fixed flavor is 

dE v — £ / dcos^P surv / dy ^Pint, (30) 

dE v J J yrn . n a uN dy 

where A is the detector's cross-sectional area, dQ Vi /dE v is the neutrino flux, and 9 Z is the 
zenith angle. The variable y is the usual inelasticity parameter in deep-inelastic scattering 
and represents the fraction of the neutrino energy going into hadronic channels, which 
means that = yE v . 

According to our previous discussion about flavor detectability, for neutral-current in- 
teractions and v T charged-currents (for energies below a few PeV), we should consider 
ymin = ythr = E^ / ' E v and y mSLX = 1, where E^ r is the threshold energy for shower de- 
tection. For v e charged-currents, y m \ a = and y m£LX = 1, since the outgoing electron also 
shower. Moreover, for v T Eq. (I30p should be modified to include the effects of regeneration. 
In fact, because of the short lifetime of tau, the v T —* r —* v T conversion takes place, with 
the result of softening the neutrino energy. In principle, one should perform a dedicated 
analysis to determine the average energy a v T has when it reaches the detector. Such a 
study is beyond the scope of our work (but see, e.g., Refs. [53,54]); we approximate this 
fact assuming that, in neutral-current events, the final neutrino energy is about one half 
of the initial energy and that, in charged-current events, the v T energy is reduced to about 
one fifth [53]. Thus, we can evaluate these integrals to obtain the total number of shower 
events: 

= £ < C +£iV s C h C (31) 

Equation (130)) is a simple form to evaluate the shower rates at neutrino telescopes. Even 
if an exhaustive discussion of the problems related to the evaluation of rates is beyond the 
scope of this work, we should mention two main sources of uncertainties. Obviously, we do 
not know the neutrino fluxes reaching the Earth's surface, neither the normalization nor 
the energy behavior. Thus, we can assume the form 

d<& v , o 



dE v 



c<S> Vi E-\ (32) 



where c is a normalization constant and <j) Ui 's are fluxes entering in the definition of the flux 
ratios in Eq. (I16p . When calculating the ratios of rates, the normalization factor cancels 
out, but the spectral index is still a crucial ingredient. In addition, the neutrino-nucleon 
cross-section is also affected by large uncertainties, mainly due to the extrapolation to very 
small x required to describe the extremely high-energy regime involved in the processes 
under discussion. 

A high-energy muon undergoes energy loss propagating in the medium, and thus, gen- 
erating showers along its track by bremsstrahlung, pair production, and photonuclear 
interactions. These showers are sources of Cherenkov light and can be easily distinguished 
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from showers previously discussed. Equation (fHUj) applies, but now 



P^{E„)ttl-e- R »f> NA <*, (33) 

where 

n^M^mt) (34) 



\a + /3Ef 

with a = 2 MeV cm 2 /g, f3 = 4.2 x 1CT 6 cm 2 /g, £? M = (1 - y) E v , and £t hr is the value 
of the muon energy to be detected. Note that the y integration in Eq. (|30|) is modified in 
such a way that y min = and y max = yf hr = 1 - Ef x jE v . 

From the previous discussion, it seems to be clear that a viable experimental quantity 
measurable at neutrino telescopes is the muon tracks N tr to shower ratio. Factorizing out 
(j) Ui from Eq. (|32p and calculating the expected rates according to Eq. (|30p . we obtain 

K = ^ = Z~ K + N^j + K ^ ■ (35) 

Note that the ratio 1Z depends on the neutrino mixing parameters through the analogous 
dependence in (f> v . . 

In order to verify these statements, we can use Eq. ([35]) once the energy thresholds 
for muons track and shower detection have been specified. To give an example, we use 
the IceCube detector setup and assume a quite conservative energy threshold E^ r = 
S* hr = 500 TeV for both types of processes in order to consider event rates well above 
the atmospheric neutrino backgrounds [55] , although for up-going neutrinos the threshold 
could be sizably smaller due to the screening effect of the Earth. For this setup, the result 
of 1Z as a function of #13 is shown in Fig. [7J in which we have evaluated ratios 1Z for the 
standard neutrino oscillation scenario and "neutrino decay I" (both hierarchies) . The error 
bars represent the statistical error associated with the number of muon tracks and shower 
events [thus, including also the uncertainties of the mixing parameters given in Eq. (|20|>] . 

For neutrino fluxes in the ratio (p Ue : <§> v : (j) Ur = 1:1:1 (standard scenario), we 
expect about 15 muon tracks against about 11 shower events in one year of data taking. 
As shown in the left plot of Fig. [7J this event rate is not enough to distinguish among 
the three scenarios considered in the text. A large overlap between the standard scenario 
and "neutrino decay I" is present for almost any value of #13. We have verified that a 
large contribution to the ratio errors is due to poor statistics. If we increase the number 
of expected events by a factor of 10 (that could be achieved in ten years of data taking), 
we observe a strong reduction of the uncertainty of the ratios 1Z, as illustrated in the right 
plot of Fig. [7J Especially, at smaller $13, the ratios corresponding to the three scenarios 
are quite separated; thus the sensitivity of IceCube seems to be enough to tell standard 
from new physics scenarios involving neutrino decay apart. Note that the effect we have 
observed increasing the event rates cannot be mimicked by a reduction of the uncertainties 
of the neutrino mixing parameters. We have calculated the ratios 7Z assuming a 10 % 
error for both #12 and 623', we have clearly observed a reduction of the errors of 1Z, but not 
sufficient to avoid the overlap among different scenarios, which means that, even for very 
precise measurements of the solar and atmospheric mixing angles, the error is dominated 
by statistics. In the opposite limit, that of infinite statistics, the right plot suggests (which 
we have numerically verified) that the distinction between the scenarios under study is 
feasible even with the current uncertainties of #12 and 623- Note that this statement has 
nothing to do with the fact that we can learn something on the neutrino mixing angles at 
neutrino telescopes, as we have extensively discussed in the previous sections. 
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Figure 7: The muon-to-shower ratio 1Z as a function of the mixing angle 9\3, evaluated for 
three of the scenarios discussed in the main text: standard neutrino framework (black boxes), 
"neutrino decay I", NH (blue diamonds), and "neutrino decay I", IH (red circles). Left plot: 
Results obtained considering one year of data taking at IceCube. Right plot: Results obtained 
considering ten years of data taking at IceCube. In both plots, the uncertainties of the neutrino 
mixing angles are those given in Eq. ft2U\) . 



5 Summary and conclusions 

The next generation of neutrino telescopes will be able to observe many high-energy events, 
thus opening the very exciting possibility to measure neutrino flux ratios. However, many 
uncertainties from both theoretical and experimental point of view, can spoil the capability 
of extracting useful information about astrophysical sources and fundamental neutrino 
properties. In this work, we have studied in detail the analytical behavior of the three 
flux ratios R e u, Rer, and R^ T , focusing on how the current errors of the neutrino mixing 
parameters can affect the theoretical expectations of R a p for any a and (3 (for flux ratios 
at the source (f>® e : <j>® : <\P Vt = 1:2:0). Due to the enormous distances traveled 
by neutrinos, the transition probabilities P a p among flavors generally assume a simple 
analytical structure. We have derived them starting from a general treatment of P a p which 
includes damping parameters to account for possible "new physics" effects in neutrino 
oscillations. We have showed that the only non- vanishing large L averaged transitions are 
those in which the non-diagonal damping factors completely disappear from the theory. 
We have then expanded the flux ratios in terms of small parameters, namely the deviation 
from maximal 623, the deviation of #13 from zero, and the deviation of #12 from its best-fit 
value. Furthermore, we have studied in detail the uncertainties of R a p connected with 
our ignorance of the fundamental neutrino physics. We have observed that the largest 
indetermination on R efl (at the level of 45 %) is due to the fact that we do not know the 
023 octant, whereas the uncertainty associated with the product cos(<5) #13 is at least a 
factor of two smaller. Due to the approximate fi — r symmetry, the same conclusion can be 
drawn for R eT . On the other hand, R^ T behaves in a very peculiar way, since the corrections 
from the standard value R^r = 1 are positive and of 0(5f 2 ). These considerations have 
been summarized in Fig. [H in which the total spread of the flux ratios caused by the 
current uncertainties of mixing parameters seems to be too large to admit the possibility 
of any realistic measurement of fundamental parameters in the neutrino sector, (e.g., 613, 
5, and the octant of #23)- 
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We have performed the same detailed study on two different scenarios involving neu- 
trino decay, considering the lightest mass eigenstate (neutrino decay I) or the lightest and 
next-to-lightest (neutrino decay II) as a stable particle(s). Once the errors of the mixing 
parameters are included, we have observed that the most promising ratio to measure devi- 
ations from the standard framework is R e ^, especially for the "neutrino decay I" scenario 
with inverted hierarchy. In that case, especially for smaller values of #13, the difference 
between the ratios can be as large as two orders of magnitude. On the other hand, the 
"neutrino decay II" scenario can be hardly distinguished from the standard result. 

Finally, we have studied the muon tracks to shower ratio 1Z at IceCube, as a physical 
observable directly connected to the experiments. We have evaluated 1Z in the standard 
picture of neutrino oscillations and the "neutrino decay I" scenario for both hierarchies. 
We have found that, contrary to the case in which 1Z is computed with the neutrino 
mixing parameters fixed to their best-fit values, the inclusion of parameter uncertainties 
reduces the possibility to distinguish "new physics" effects in neutrino oscillations. Thus, 
we believe that any realistic analysis of physics reach of neutrino telescopes concerning flux 
ratios should carefully include parameters uncertainties. 
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A Series expansion formulas 

A.l First-order series expansions of neutrino flux ratios with 
two non-zero normalized damping factors 

The neutrino flux ratios using two of the normalized damping parameters equal to 1, i.e., 
d\ = d,2 = 1, and the last one equal to zero: 

(36) 
(37) 
(38) 

The neutrino flux ratios using two of the normalized damping parameters equal to 1, i.e., 
d\ = ds = 1, and the last one equal to zero: 

o / 7 , n 4cos 2 (M 16sin(2g 12 ) 32cos 2 (g 12 ) 

R e ^ dl =d 3 = l) = 3 _ cos{2eu) ~ [3 _ cos(2M] 2*i2 " [3 - CO s(2M]^ 23 ' { ) 

t> (j j u 4cos 2 (fl 12 ) 16sin(2fl 12 ) 

K eT {di — d 3 — 1J — — — - - j- — — 775012 

3 - cos(26q 2 ) [3 — cos(26 l i 2 )J^ 

32cos((5) cos 3 (9 12 )sm(6 12 ) _ 8sin 2 (2fl 12 ) 
[3-cos(2# 12 )] 2 [3 - cos(2# 12 )] 2 



Re/i(di 


= d 2 


= 1) 


= 2 + 2 cos(<5) sin(40 12 )#i 3 + 2[3 + cos(40i 2 )]<5 23 


Rer{dl 


= d 2 


= 1) 


= 2 + 2 cos(5) sin(40 12 )# 13 - 4sin 2 (20 12 )5 23 , 


Rfir(di 


= d 2 


= 1) 


= 1-4(523- 
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o ( , , u 1 , 4co S (^) S in(2g 12 ) fl 8cos 2 (M 

i? MT (di = d 3 = 1) = 1 + —5 , . gi 3 + o T^-^rfe- 41 

3 — cos(2t^i2j 3 — cos(2#i2j 
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